Opposite effect of spin-orbit coupling on condensation and superfluidity 
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We investigated effects of a Rashba-type spin-orbit coupling (SOC) on the condensed density and superfluid 
density tensor of a two-component Fermi gas in the BCS-BEC crossover at zero temperature. In anisotropic 
three dimensions (3D), we found that SOC has an opposite effect on condensation (enhanced) and superfluidity 
(suppressed in the SOC direction) and this effect becomes most pronounced for very weak interactions and the 
SOC strength being larger than a characteristic value. Furthermore, as functions of SOC strength, the condensed 
density changes monotonously for all interaction parameters while the superfluid density has a minimum when 
the interaction parameter is below a critical value. We also discussed the isotropic two dimensional (2D) case 
where analytical expressions for the gap and number equations were obtained and the same phenomena was 
found as that of the 3D case. 
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Introduction. Spin-orbit coupling (SOC) is a central topic 
in condensed matter physics. First, it plays an essential role 
for the realization of nontrivial topological states which are 
discussed intensively nowadays IjJ]. Secondly, as was shown 
by Gor'kov and Rashba f^, SOC can induce a nontrivial 
spin-triplet pairing field which leads to significant changes 
in the properties of superconductors |3]. Quite recently, ef- 
fective SOC was realized for bosonic ultra-cold atoms 
by dressing two atomic spin states with a pair of lasers |4]. 
With the anticipation that this novel technique is also applica- 
ble to Fermi atoms, practical proposal of generating SOC in 
fermionic '^'^K atoms with tunable interaction through Fesh- 
bach Resonance is given in (sl]. 

Motivated by this new progress, effect of SOC on the pair- 
ing and superfluid nature of Fermi systems in the BCS-BEC 
crossover has become a cutting-edge field recently because 
of its broad interests in condensed matter physics. The spin- 
triplet pairing fields and anisotropic nature of the superfluid- 
ity induced by SOC were investigated in |6] and proposal for 
detecting this phenomenon was given in |7] through measure- 
ment of the momentum distribution and single-particle spec- 
tral function. On the other hand, SOC significantly enhances 
the pairing phenomena as was shown by the exact two-body 
solutions 1 8] where a new bound state (rashbons) emerges and 
many -body mean-field calculations llo Ul ill . 

In this Letter, we study the effects of SOC on two funda- 
mental quantities: condensation and superfluidity. Conden- 
sation is well described by the concept of off-diagonal-long- 
range-order |fT3]. However, Landau's approach of calculation 
of superfluid density (tensor) is only applicable to systems sat- 
isfying Galilean transformation 11311 . For systems in the pres- 
ence of SOC obviously violating Galilean transformation, we 
gave the general method of calculating the superfluid density 
tensor Furthermore, we found that at zero temperature, SOC 
enhances condensation while suppresses superfluidity in both 
3D and 2D. Up to our knowledge, this is the first demonstra- 
tion of such opposite behaviors of condensation and superflu- 
idity driven by SOC and renews our previous knowledge that 



these two phenomena change in the same direction with other 
influencing factors (such as temperature and disorder). 

The model.ln the presence of SOC, the system 
of a two-component Fermi gas can be described 
by the finite temperature grand-partition function 
Z = /d[?/5^,V'CT]exp(-5[?/^CT,i/'cr])_ (ft = fc_B = 1 
through out) where the action S[ipa,tpa] is given by 
S[^p,,^^] = Jl^dT f d^rJ2J^„drJja +Ho + Hi] with 
/3 = 1/T , (7 =^;\r denoting spin, 4'cri4'cr being the Grass- 
mann fields, and d (= 2, 3) being the dimension. We focus on 
a Rashba-type SOC [2*] and the single particle Hamiltonian 

density can be written as ijj)—'tp (j,p + T-Lsoj V' where 
ill — [V't, V'J.]"'" is the collective fermionic field, the kinetic 
operator = p^/(2to) — with ^ being the chemical 
potential, the Rashba term "Hso = A (ct x p)^ with a being 
the Pauli matrices and A being the SOC strength. The 
singlet-channel attractive interaction can be characterized by 
a contact interaction parameter g (< 0) and correspondingly 
Hi = gip^Tp^Tp^tp^. 

In order to study the Fermi-pairing problems, we em- 
ploy the Hubbard-Stratonovich transformation [14] to can- 
cel the four-body interaction term Hi by introducing a pair- 
ing field A (r, r). After integrating out the fermionic fields, 
we obtain the effective pairing action as S'e// ['^i^] — 
- /(f dr J d'^r I A (r, T)f /g ~ l/2Tr In {Q^^\ where the in- 
verse propagator t/^^ is 



dr+^p % A 

7^ dr+ip -A 

-A dr-ip % 

A 7p dr~ip 



(1) 



with7p = \iPy + ipx)- 

Mean-field theory. At the mean-field level A (r, r) = 
Ao which is referred to as the gap parameter, the effec- 
tive pairing action becomes Seff [A, A] = —pVAl/g — 
V2Ep.^^„ In [det^p'^LJ where G'}^^ is the momentum- 
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frequency representation of Eq. ([T]i with A (r, r) = Aq, V is 
the size of the system and tj„ = (2n + 1) 7r//3 are the Fermi 



Matsubara frequencies. From det Q 



0, the excitation 



spectrum can be obtained as Ep^s = \J (Cp + <^ l7p 



and -Ep ^ — —Ep^s where — €p ~ fi with Cp = p^/2rn 
and 5 ~ ±1 is called helicity. Finally, by using the thermo- 
dynamic relation ft — — l//31nZo, we have the thermody- 
namic potential fin = -VA^/g + l/2Y,pj (fp - Ep^s) - 
l//3^p hi (l + e^'^^p "') from which the gap and num- 
ber equations are given by 
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Eq. ^ and Eq. ^ are the generalized BCS gap and num- 
ber equations in the presence of a Rashba-type SOC which 
have been investigated in detail to study the ground state and 
finite temperature properties of this novel system. The key 
discovery is that the increased density of states by SOC plays 
a crucial role for the understanding of the pairing enhancing 
phenomena 11 ill . With these results in mind, we now move 
on to the calculation and discussion of condensed density and 
superfluid density tensor 

Condensed density.For the Fermi pairing problems, 
the condensed density is generally defined as lll2ll 

Tic = l/"^I]p,ss' |('0p,sV'-p,s')l ■ *^ system 

considered in this paper, the singlet-channel attrac- 
tive interaction supports a singlet-pairing field while 
SOC can simultaneously induce a triplet compo- 
nent. Within the mean-field theory, spin-singlet and 
-triplet pairing fields are given by iQ]: (V'p.t^'-p.i) = 
AoE5tanh(/?£;p.5/2)/(4i;p,5) and (V^p.tV'-p.t) = 
-Ao (7p/ ITpI) '5tanh {l3Ep^s/2) / i4.Ep^s), respec- 
tively. The spin-singlet contribution to the condensed fraction 
was first discussed in |15] where it was shown to behave 
non-monotonously with a minimum as a function of SOC 
strength for weak enough interaction parameter. In this Letter, 
we take both pairing components into consideration and the 
full condensed density becomes 



A2 
^0 



tanh" 



(^) 



(4) 



At zero temperature, repulsive interactions between Fermi 
pairs (Bosons) result in depletion of the condensate which is 
a familiar phenomenon for interacting EEC systems. 

Superfluid t/eni/fy. Unlike the condensate density, the super- 
fluidity is a kinetic property of the system. By Landau's theory 
11131], the normal mass of the system can be obtained through 



the calculation of the total momentum carried by excitations 
when the system is enforced in a uniform superfluid flow with 
velocity 



p = ^p/ (ii;p,, - p • V J , 



(5) 



where cr is a conserved quantum number which is spin in the 
absence of SOC, f {x) — 1/ [e^^ ± l) is the Fermi/Bose dis- 
tribution function depending on the nature of the excitations, 
and Ep^rr ^ P • is the excitation spectrum for moving sys- 
tem obtained from Galilean transformation. At zero tempera- 
ture, no excitations are created at very small and the super- 
fluid density coincides with the total density. 

However the situation is dramatically changed in the pres- 
ence of SOC where Galilean transformation is violated. In 
order to calculate the response of the system to a uni- 
form superfluid flow in the presence of SOC, instead of us- 
ing Eq. (|5]l which is no longer valid, we calculate the 
increasing of thermodynamic potential Sfl = il (vg ) — 

Q,{0) = {l/2)Vj2ar]'"T''"'S,ariVaVn from which ns,arj = 

1/ (mV) [9^r2 (vs) /dvs,advs,,-i]^ _^ is defined as the super- 
fluid density tensor lfl6ll . A convenient approach of generat- 
ing such superfluid flow is applying a "phase twist" to the or- 
der parameter ||l7ll : A (r, r) = Aoe"i and correspondingly 
Vs — q/2m. Therefore the superfluid density tensor can be 
defined as Us^ct] = 4™ [d^il (q) /dqadq,^] 

By substitution of A(r,r) = Aoe'i"" into Eq. ([T|, 
the thermodynamic potential for moving system can be ob- 
tained as n (q) = ~VAl/g + Ek (^k - k • q/2TO^ - 

g/3(Bp.,-k-q/2m)" 
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l/(2/3)Ep,=i.4ln 

+ q^/ (8m) and Epj are solutions of 



Ei 



£^p,»7q/2 - ?p7p 



where ^p 



^ 0, (6) 



withwp = f^ + A^ + |7p| -|7q/2| , ana Ap = Jm (^7p7q^ 
In the presence of SOC, when the whole system is moving 
with a uniform velocity, the original four excitation spectrums 
(Ep.s, —Ep s) are strongly coupled and correspondingly Eq. 
(|5]l is not well defined now. Superfluidity in systems that vi- 
olates the Galilean transformation has also been discussed in 
the bosonic systems in the presence of SOC where the criti- 
cal velocity has been discussed with the same method used in 
the calculation of the excitation spectrum for moving systems 
ill. 

Combined with Eq. calculation of the second-order 
derivative of ft (q) with respect to qi is straightforward al- 
though tedious. Finally, the superfluid density tensor is ob- 
tained as 



\2m/ 



(7) 
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S^p hp\ Ep^s 



where Yp^s = Pf (Ep^s) [1 - / (Ep^s)] with / (x) being the 
Fermi distribution function, Ug^xx — ng^yy — ||, and 
Us^a^jj = 0. In 3D, the anisotropic nature of the superfluid 
can be evidently seen from Ug.zz 7^ Since SOC does not 
affect motion in z direction, spin is a conserved quantum num- 
ber and Us.zz has the same form as that obtained from Eq. (|5]l. 
However, superfluid motion in x, y direction is dramatically 
changed by SOC. Most interestingly, a new term (last line in 
Eq. dHJ) emerges due to SOC and it is not zero at T = which 
means suppression of superfluidity in x, y direction. The first 
line of Eq. (|8]l can be understood in spirit of Landau's theory 
where momentum carried by excitations Ep s is now shifted 
by SXm due to SOC. At T = Tc where Aq = 0, both Us^z 
and II equal to zero which is crucial for the correctness of 
our results. In 2D, the system is isotropic where the superfluid 
density is only given by ||. 

With the anticipation that the mean-field theory is sufficient 
to capture the qualitatively correct physics in the whole BCS- 
BEC region at zero temperature as demonstrated in [19] in 
the absence of SOC, we only focus on the zero temperature 
behaviors of the condensed and superfluid density. 

Results and discussion.Because, at zero temperature, be- 
haviors of the condensed and superfluid density (risn) are the 
same in 3D and 2D cases with only quantitative differences, 
we will discuss the 3D case in detail and give a brief discus- 
sion on the 2D results. 

As usual we regularize the contact interaction parameter g 
in Eq. ^ by the experimentally related scattering length a 
through 1/g = m/ (Ana) -l/VJ2p^/ (2ep)- With the gap 
and chemical potential obtained from the self-consistent so- 
lutions of Eq. ^ and Eq. (O, we numerically calculate Eq. 
dUl and Eq. (|8]l and the results are shown in Fig.l. As can 
be seen from Fig. [Tta) and (c), the condensed density is al- 
ways enhanced by SOC. Nevertheless, seen from Fig. ^a), 
we can still define a characteristic value roughly located at 
Ac ~ 0.5vf where vp — kp/m and kp is defined through 
kp — Stt^u. For A < Ac, condensed fraction defined as 
fic/n, increases only slightly which can also be seen from Fig. 
[Tfc) that the solid orange and dotted blue lines almost coincide 
with the dash black line (where A = 0). Only when A > Ac 
can ric/n have a significant increase. For 1/kpa — !> +00, 
the effect of SOC becomes very weak and nc/n — > 1. For 
Ac ^ vp, we also have ric/n — 1 which agrees with the 
statement that SOC can produce a bound state and thus in- 
duces a crossover even for very weak interactions lIsl lTlll . 

On the contrary, the superfluidity is always suppressed by 
SOC as can be seen in Fig. [TJb). Furthermore, as a func- 
tion of A, superfluid fraction || /n varies non-monotonously 
with a minimum located at A'c (denoted by red points in Fig. 
[Tfb)). This minimum exists for interaction parameter below 
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FIG. 1: (Color online) Condensed and superfluid fraction as func- 
tions of \/vf and l/ikpa) in 3D. In (a) and (b), l/(fcFa) are given 
as: —2, —1.2, —0.4 and 0.4 for the same line species from below. In 
(c) and (d), \/vf are set to be 0.4 for solid orange lines, 0.6 for dot- 
ted blue lines, 1.6 for dot-dashed red lines and 2.0 for media-dashed 
green lines. Dashed black line in (c) corresponds to A = 0. Red 
points in (b) and (d) are explained in text. 



a well-defined critical value l/kpUc « —0.13 which can be 
determined in Fig. [TJd) as the rightmost red crossing point. 
However, we find that SOC never destroys the superfluid com- 
pletely for all interaction parameters. The minimum super- 
fluid fraction (^s,||/"■)J-f^j„ ~^ 0.5 for 1/kpa — > ~oo. When 
1/kpa > 1/kpac, this monotonic behavior disappears and 
||/n decreases with increasing A. 

The opposite behavior of condensation and superfluidity 
controlled by SOC strength is found to be most pronounced 
for very weak interaction parameters and A > max(Ac, A^). 
As can be seen from the solid orange lines in Fig. [Ha) and 
(b), for A — 2vp, the condensed fraction ric/n is increased 
by 0.78 and the superfluid density is suppressed by 0.38. For 
A < max(Ac, Ac), the superfluid density decreases quickly to 
its minimum value while the condensed density changes only 
very slightly. 

In 2D, divergence of the integral over momenta can be 
cured by substituting 1/g = -l/V X)p V (2ep + 
into Eq. (|2]i where Eb is the binding energy and be- 
comes the controlling parameter for the BCS-BEC 
crossover problem. With dimensionless parameter 
given by A = A/uf, A = 
Eb = Eb/Ep where Ep = 
fined through kp = 2im, Eq. 



fi/Ep, Ao = Ao/ Ep, 
kp/2m and kp is de- 
^ and Eq. ^ have 



analytical expressions as 1 = 



/2 



X' 1 



H[fi,Ao,X 



InEb 
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While the nonzero normal density at zero temperature induced 
by SOC is a direct result of the triplet pairing field. In |9], it 
is demonstrated that, at zero temperature, this triplet pairing 
field induces a non-zero spin susceptibility which implies a 
residual normal fluid 1.2 3il . An interesting future work is the 

combined effect of disorder and SOC on these two phenom- 

1 — 

ena where a bose-glass state [21] may show up based on the 
conclusion that disorder and SOC both suppress superfluidity 
while depletion of the condensate induced by disorder may be 
weakened/cancelled by SOC. 
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Note Added: After finishing this paper, we note that the 
full condensed density is also discussed in ll24ll . 



FIG. 2: (Color online) Condensed and superfluid fraction of 2D 
Fermi gas as functions of A/uf and Eb / Ef- In (a) and (b), Eb /Ef 
are given as: 0.001 for solid orange lines, 0.01 for dotted blue lines, 
0.1 for dot-dashed red lines and 1.1 for medium-dashed green lines. 
In (c) and (d), \/vf are set to be 0.4, 0.8, 1.2 and 2.0 for solid or- 
ange, dotted blue, dot-dashed red, medium-dashed green lines, re- 
spectively. Dashed black line in (c) corresponds to A = 0. Red 
points in (b) and (d) are explained in text. 



= In 



A2 

^0 



/i 



ii'l/i, AqjA) respectively 



with H yjl, Aq, Aj and K yji, Aq, Aj given in ||20|]. Results 
of superfluid and condensed fraction are shown in Fig. 2. The 
same phenomena discussed in 3D are also found in 2D case. 
However, we find that Ac drifts leftwards when increasing 
Eb I Ep as can be seen from Fig. I2a). The critical interaction 
parameter for the appearance of a minimum point of nsj|/n 
is Eb « OMSEf- 

Conclusions. \n summary, general formulae were obtained 
for the condensed density and superfluid density tensor of a 
two-component Fermi gases in the presence of a Rashba-type 
SOC. At zero temperature, we found that superfluidity in the 
SOC direction is suppressed while condensation is enhanced 
by SOC and this phenomenon becomes most pronounced for 
very weak interaction parameters and A > max(Ac, A^). Fur- 
thermore, the superfluid fraction exhibits a non-monotonic be- 
havior with a minimum as a function of SOC strength when 
interaction parameters is below a critical value while the con- 
densed fraction increases only monotonously with either inter- 
action parameter or SOC strength. These phenomena happen 
in both 3D and 2D cases. 

Finally, we point out that there is an essential difference 
considering the mechanism of suppressing superfluidity by 
disorder 021 



2211 and SOC. Superfluid motion is suppressed 



through energy dissipation due to scattering with impurities. 



[1] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005); 

M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010); 

X. L. Qi and S. C. Zhang, ar Xiv: 1 108.2026 T Rev. Mod. Phys. 

(to be published)]. 
[2] Y. A. Bychkov and E. I. Rashba, J. Phys. C 17, 6039 (1984); 

L. P Gor'kov and E. I. Rashba, Phys. Rev. Lett. 87, 037004 

(2001). 

[3] F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Rev. Mod. Phys. 

77, 1321 (2005). 
[4] Y.-I. Lin et at.. Nature (London) 471, 83 (201 1). 
[5] I. D. Sau etal., Phys. Rev. B 83, 140510 (2011). 
[6] I. P. Vyasanakere, S. Zhang, and V. B. Shenoy, Phys. Rev. B 84, 

014512 (2011) 

[7] Hui Hu et al. , larXiv: 1 105:2488] IPhys. Rev. Lett, (to be pub- 
Ushed)]. 

[8] 1. P Vyasanakere and V. B. Shenoy, Phys. Rev. B 83, 094515 
(2011). 

[9] Li Han and C. A. R. Sa de Melo,larXiv:1106.3613J 
[10] M. Iskin and A. L. Subasi, Phys. Rev. Lett. 107, 050402 (201 1). 
[11] Z.-Q. YuandH. Zhai, arXiv: 1105.2250 

[12] C. N. Yang, Rev. Mod. Phys. 34, 694 (1962); A. I. Leggett, 
Quantum Liquids. Base Condensation and Cooper Pairing in 
Condensed-Matter Systems (Oxford University Press, Oxford, 
2006). 

[13] E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Part 2 

(Butterworth-Heinemann, Oxford, 2002). 
[14] A. Altland and B. Simons, Condensed Matter Field Theory 

(Cambridge University Press, Cambridge, England, 2006). 
[15] L. Deir Anna et al. arXiv: 1 108. 1 132 i'l. 

[16] E. Taylor, A. Griffin, N. Fukushima, and Y. Ohashi, Phys. Rev. 

A 74, 063626 (2006); N. Fukushima, Y. Ohashi, E. Taylor, and 

A. Griffin, ibid 75, 033609 (2007); Edward Taylor, Phys. Rev. 

B 77, 144521 (2008). 
[17] M. E. Fisher, M. N. Barber, and D. lasnow, Phys. Rev. A 8, 

nil (1973). 

[18] Qizhong Zhu, Chuanwei Zhang, and Biao Wu, 
arXiv: 1109.58 11 

[19] G. E. Astrakharchik, I. Boronat, I. Casulleras, and S. Giorgini, 
Phys. Rev. Lett. 95, 230405 (2005). 



5 



[20] H(^fl,A,X] = \A'^[N{29,K.)-M{2d,K)]/ /jf l/\/l - Ksin^ xdx are the first and second kind in- 

complete elliptic integrals respectively with < <^ < tt in this 



1/4 



, K[fi,A,X) = \F{29,k)/ paper. 

[21] K. Huang and H. F. Meng, Phys. Rev. Lett. 69, 644 (1992). 



(xl + aA where xo = A + A^ 2k = 1 + xo/Jx^ + A'^, [22] G. Orso, Phys. Rev. Lett. 99, 250402(2007). 

{ , - \i/4 [23] A.J. Leggett, Phys. Rev. 140, A 1869 (1965). 

tanS = 'X/lxl+A"^] ,N{(I>,k) = E{(j),K) + {K~ [24] L. Dell' Anna e/ a/., arXiv: 1108. 1132^-3; Beibing Huang, Shao- 



1)F(0,k), and M(0, k) ^ Ksin ^cos0/v/l-Ksin^0. long Wan,| arXiv: 1109.39701 

-B (<^, k) = /jj** \/ 1 — Ksin'^ xdx and F ((/>, k) — 



